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ABSTRACT 


шаны Thesis computerizes aircraft Naval Air Training and 
Me@etrating Procedures Standardization Program (NATOPS) flight 
performance charts and shows it to be feasible with a high 
degree of accuracy. The computer programs developed are 
ШОО lake-ort and Cruise Performance for the A-6E aircraft 
and are adaptable to hand-held calculators, desk calculators, 
n existing aircraft systems computers. The feasibility, 

ШОШО тес, апа techniques shown in this thesis are applicable 
понапу aircraft, both fixed and rotary wing. 

"much improved and more accurate method of mission plan- 
SS heeded to reduce NATOPS chart errors during pre-flight 
mins and to reduce the accident potential of air crews who 
make performance decisions during flight from past experience 
Са intormation, without accounting for all variables which 
Mereet aircratt performance. A large percentage of Navy/ 
Marıme Corps aircraft accidents has been directly or indirectly 
led to misinterpretation of NATOPS charts and/or poor 
нениноштапсе decisions. Computerized NATOPS performance charts 
fewest Oniticantly reduce human errors inherent to visual 
chart interpretation and, when coupled with aircrews having 
ШОШ r access to make airborne aircraft performance deci- 
= based on all performance variables, will greatly enhance 
safety of flight in completing any or all mission phases, 
thereby increasing operational readiness of all Navy/Marine 


orps fleet units. 
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шас а аса а роі 


objective of this thesis was to investigate the feas- 
MONEY Ot computerizing aircraft Naval Air Training and 
Seerating Procedures Standardization Program (NATOPS) flight 
Merrormance charts with a high degree of accuracy and to 
develop computer programs which are adaptable to hand-held 
mmpemretors, desk calculators, or existing internal aircraft 
Stems computers. The NATOPS performance charts of opera- 
aircraft used in the Navy and Marine Corps are a 
conglomeration of graphs, charts, curves, and tables which 
Ire quentiy cumbersome and very susceptible to user error 
meen: Interpretation and manipulation. The charts present 
p ooPmance information for all aspects of flight operations, 
ENS ke-off to final landing. Therefore, it is imperative 
ent crews be able to interpret and utilize these 
mers tO obtain optimum and efficient performance from their 
тесту е aircraft in all phases of the specified mission. 

Mission planning is an essential and key element of flight 
Operations and must insure a high degree of safety and maximum 
Mt zation of the aircraft performance in order to success- 
my accomplish the mission. This planning requires use of 
the NATOPS charts to determine take-off distances, climb 
profiles, cruise distances, low-level distances and fuel used 
an all mission phases. Such mission planning is very time 


Em ло отш ету susceptible to error since it requires 








KJ] interpolation between curves which are as wide as 

the user's pencil lead, and various information is transferred 
ИП спатс to chart. For example, even with extremely precise 
нини зегиі зпсегргебсасіоп of take-off distance, an error of 
IN Eeet could be made due to pencil width alone, exclusive 
of interpolative errors. The combination of these human and 
mechanical errors could build to as much as 500 feet, which 
would be very critical when taking off at maximum gross weight 
on a hot day. Errors in fuel management would be equally as 
OI if the mission profile dictated return to base or 
ПСВ minimal fuel. Also, it is noted that utilization 
etme NATOPS performance charts is generally limited to pre- 
flight planning and not feasible for use in high performance 
Penerart, and in several of the helicopters, for in-flight 

Pm es in mission profile. A NATOPS pocket check list is 
N TO! In-flight planning, which is small and therefore 
un C to even higher human and mechanical interpolation 
етого. 

Therefore, a much improved and more accurate method of 
Mission planning, both before flight and during flight, is 
ЕЕ A desk top calculator or hand-held computer with 
pre-programmed chips could eliminate these errors during 
EL "vsht planning, while in-flight errors could be reduced 
substantially with performance equations programmed into an 
on-board systems computer of the aircraft, or a hand-held 
Eomputer with a dedicated or read only memory (ROM), a pro- 


mem chip, Capable of storing performance equations for all 








mission phases. Not only would computerization reduce the 
human errors of flight crews who have the NATOPS charts avail- 
EN Ut would reduce the accident potential of flight crews 
who are airborne without NATOPS and make performance decisions 
ШЕШ рас: experience. Such decisions, without accounting for 
Eu ariables which affect aircraft performance, could lead 

EN c: errors in predicting the capability of the aircraft 
KS ely complete a mission phase; and have been cited as 

the Cause of a large percentage of Navy/Marine Corps accidents 
Ma tixed and rotary wing aircraft. Computerization would 
also Ey enhance shipboard operations by providing in- 
Bene Bingo (minimum fuel) profiles to both flight crews and 
the air base, reducing fuel management errors based on exper- 
Mence and interpolation of charts. 

ШЕ thesis computerizes A-6E aircraft NATOPS take-off and 
Eu performance charts [Ref. 1]. However, the feasibility, 
procedures and techniques are applicable to any aircraft, 
both fixed and rotary wing. The feasibility aspect of the 
EL Ss € continuing application of thesis research con- 
meted гог Computerization of A-7E performance [Ref. 2]. The 
EN cams are written in BASIC language for the Hewlett-Packard 
moun computer, but knowledge of the language is not required 
EN cuonstrate the feasibility of accurate performance com- 
puterrzation. The performance equations within the programs 
sS directly utilized In most scientific computer languages, 
With easy adaptation to any computer such as hand-held, desk 


mop, Or aircraft systems computers. 
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SES PERIMENTAL PROCEDURE 


A. GENERAL ORGANIZATION 

Two separate, but equally important, A-6E aircraft per- 
formance charts were evaluated for computerization feasibility. 
Mier tarst chart, "Take-off Performance" [Ref. 1, Fig. 11-11], 
Was Utilized since its family of relatively simple, or low 
order polynomial equations, provided use of basic curves in 
Studying and developing the programming techniques. The 
EE Chart, "Cruise Performance," [Ref. 1, Fig. 11-110], 
provided a more complex familv of curves which required more 
extensive analysis in program development. Empirical data 
n K on from the specific NATOPS chart and the Least-Squares 
method of polynomial approximation was used to determine the 
coefficients of each curve in the respective family of curves. 
A polynomial expression was then approximated using the same 
Ш ПОП tO determine the family of coefficients. In other 
Words, polynomial expressions were determined which approxi- 
mated the coefficients of any curve which would fall within 
ЕПЕ bounds of the respective family of curves. This technique 
precluded the need for interpolation between curves. A 
detailed example of the Least Squares method is presented in 
Appendix A, while Figure 1 presents a simplified example of 


approximating a family of curves. 


el 





101) 





XC 1) 


Басите І. Family of curves. 


Teast Squares approximations yield polynomial expressions 


E each curve in the form of 


Р 2 
EN S ао ХІ 
Y(2) = a + a аа: X 

20 1 2 Pr 2 
Y(3) = a — Хеба х,“ 
30 spi 5 $2 3 
Y(4) =a + a хх зна a 
40 41 4 42 4 


where for a.,. 
1) 


"Xu, iS The Specific curve and j is the 
Moetficient subscript. 
Pwepolynomial expression for each coefficient family is 
determined in the form of: 
ACOD = pu v by Sue b, e C 
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mere "UC" is the curve variable. 
mens the expression tor the curve passing through point 


En Figure 1 for X; = С, с Шен 


5) - 00) + A() * X, « * AC) * X, „2. 


24519 

The Least Squares method was programmed for the Hewlett- 
Packard 9850A computer complemented with the Hewlett-Packard 
9266A plotter which was used to verify that the plot of the 
polynomial expressions did, in fact, duplicate the specific 
WETOPS curve. The accuracy of the polynomial approximations 
was determined by evaluating the percentage of deviation of 
the approximated curve at specific points from the same points 
mm the actual NATOPS curve. It was possible to achieve ваше аве 
than ninety-seven percent (97%) overall accuracy. 

Accuracy of the polynomial expressions for each coefficient 
family was determined by comparing the computed value of the 
coefficient evaluated at each curve of the family with the 
@ereinal coefficient approximated for each curve. А more 
oiled description of the computer techniques is contained 


Mma Ref. [2]. 


B. TAKE-OFF PERFORMANCE 
The A-6E aircraft NATOPS normal take-off performance chart 
(Fig. 3] was photographed and enlarged five times its original 


Size to allow more accurate measurement of empirical data. 
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The chart was then separated into sub-charts composed of 
Families of curves for individual evaluation of computer- 
ization feasibility. The dependent variable or output of 
KS first sub-chart became the independent variable or input 


NE succecdaing sub-chart. Table I shows the relationship 


ENUsub-charts to the independent and dependent variables. 


TABLE I 


шара UOTT PERFORMANCE SUB-CHARTS 


Віа с е Altitude Temperature 


Gross Weight Reference 
Index Take-off 
Distance 
Take-off Distance Wind Take-off 
Distance 
Empirical data were taken from the respective sub-chart, 
the coefficient of each curve was determined, and the family 
E Wives was plotted to verify accuracy with the original 
E liri. Once the points verified for each curve were 
DT One percent of the original empirical data, the coef- 
ficients of the family of curves were determined. A computer 
program was written and the coefficients were verified for 
icy within two and one-half (2-1/2) to three (3) percent 
Or the coefficient previously obtained for each curve. The 


independent variable for the coefficients of the family of 


Bes is the specific "family of curves parameter," such as 


14 





Pressure altitude in Sub-Chart I. The dependent variables 
шаце coefficients of the curve passing through the specific 
maby OL curves parameter. These coefficients are then used 
Не output equation of the specific sub-chart. 

i ESub-Chart I, Pressure I pide iE was relatively 
E tO determine the coefficients of the family of curves, 
Since there were the same linear relationships between each 
Eu But Sub-Charts II апа III, "Gross Weight" and "ТаКе- 
off Distance with Wind," had to be broken down into three and 
two sub-families of curves, respectively. This was required 
вашае the non-linear relationship between curves would pro- 
Sueceimaccurate coefficients for the family of curves. Ву 
Overlapping smaller sub-families of curves whose relationships 
mere more linear, the desired accuracy could be achieved. The 
range covered by each sub-family of curves was determined by 
 Ріа! апа еггог iteration until the desired accuracy was 
al c. Coefficients for the sub-family of curves would 
Ku sed in the output equation for the specific sub-chart. 
SD amily coefficients used depend upon the range in 
The desired family of curve parameter would lie. The 
шашшысет program was written to check for the limits of each 
sub-family range. 

E|ueccsatystactory programs were written for each sub-chart, 
they were integrated into one computer program which would 
Erud a final output of take-off distance with wind. The 
final program also includes equations which calculate take-off 


Deity, ahd unsafe and not-recommended take-off distances. 


185 





The coefficients for these equations were determined in the 
Same manner as individual curves for the sub-charts. A wind 
component computation using basic trigonometry identies was 
ШЕР (0 input winds into Sub-Chart III, and a round-off 
routine was used prior to the final output printout. The 
Emil orff criterion established for take-off velocity was 
one-half (0.5) knots, while that for take-off distance was 
twenty-five (25) feet. Twenty-five feet was used for take- 
Off distance round-off, since at high take-off velocities 
this value would be negligible in rounding down, but would 
still provide агас от ОЕ safety tor round-up to the nearest 


one hundred (100) feet. 


Mee CRUISE PERFORMANCE 

The A-6E NATOPS Cruise Performance Charts are a series of 
Separate charts or four phases which require that the output 
Ot each chart be transferred as the input to the succeeding 
a until arriving at the final output of "fuel flow." 


ese four phases are: 


Phase 1 Clean Aircraft Transfer Scale 
Phase II Aircraft Reference Number 
Phase III Pounds of Fuel per Nautical Mile 


Phase IV Fuel Flow 


Since Phase I chart consisted of the most non-linear 
Eves, in contrast to the linear curves of Phases II, III 
IN, and in view of project time constraints, only Phase I 


was evaluated for computerization feasibility. This chart 
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К 150 photographed and enlarged to facilitate extraction 
ШОО rate empirical data. As in take-off performance, the 
Enc was separated into sub-charts for individual evaluation 
ЕЕЕ ШІ 1165 of curves. This technique varied from that of 
take-off performance in that Sub-Chart II had an overlying 
EU x Curves which had to be treated as separate entities. 
The dependent variable or output of one curve family became 
an independent variable or input to the second curve family. 
ШЕГЕ ТТІ presents the sub-charts with respective curve 


тез ап variables. 


TABLE II 


ST TT CRUISE PERFORMANCE SUB-CHARTS 


mS Chart | Family of Curves Parameter Input Output 
Т Pressure Altitude Gross Weight | Mach 
TET cele Mach Шаси 
Scale 
IIT ла Count Træn ter Mach 


Scale 





termination of the coefficients of each curve, accuracy 
Mrittcation, and determination of curve family coefficients 
were accomplished using the same method described in Take-off 
ESS mance. Sub-Chart I, "Pressure Altitude," and Sub-Chart 
II, "Guideline," had to be broken down into two overlapping 
sub-families of curves to achieve the desired accuracy. The 
range of each sub-family was determined by trial and error 


Until the desired accuracy was achieved. Coefficients for 


П 





Saemeranmily of curves were used in the output equation for 
E puter program written for each family of curves. Each 
Bam with sub-families was checked for the range in which 
the family of curves parameter would lie to ensure that the 
Meen Coefficients were utilized. 

Smee Sub-Chart 11 had overlying curve families, which 
Ned that the input variable, Mach number, would intercept 
шопідїеїїпе at the zero drag count, and proceed down that line 
untill intersecting the desired drag count line, a comparing 
and decrementing routine had to be developed. A graphical 


example of this routine is shown in Figure 2. 


DRAE 
(004T LINES 
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Figure 2. Mach comparison and decremation. 
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тивно routine, initial Mach (MØ) from Sub-Chart I is 
miput, as the ee pendenesvarsable, into the zero drag count 
Ex on which outputs the dependent variable, initial transfer 
Beale (TP). Tf then becomes the independent variable which 
determines the specific guideline to be used by establishing 
Me ricients. The independent variable іп the guideline 
Eu en is Mach, which calculates the dependent variable, 
Bester Scale (TI). Tl is then input as the independent var- 
u III O the selected drag count equation which calculates 
miemaependent variable, Mach (Ml). МІ 15 then compared to 
L S u Ous Mach which was input to the guideline equation. 
For the first iteration, Ml is compared with Mf. If the 
ШҮСІПГЕ value of the difference between the two Machs is not 
Hess than or equal to 0.01, then Mf is decremented, А ғ .01, 
Па шем transfer scale, Tl', is calculated and input back 
Ehe drag equation to output's new Mach, Ml'. This iter- 
ative and comparison routine is continued until the comparison 
parameter is satisfied. The compared Machs are then averaged 
and output as Best Cruise Mach, M2, while the last computed 
Bs output as Transfer Scale, T. M2 and T are the final. 
puis of Sub-Chart II. 

mee Satisfactory programs were written for each sub-chart 
Beemer i tied for accuracy, they were integrated into one 
puter program to yield the final output of Best Cruise Mach 


and Transfer Scale. The numerical values of these outputs 


ere not rounded off for aircrew use, since they would be 


nS 





transferred as inputs to subsequent phase charts of the cruise 


performance computer program in follow-on thesis work. 
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ЕГІ БЕЗШЕІТ5 АМ) DISCUSSION 


fee TAKE-OFF PERFORMANCE 
MSE NATOPS normal take-off performance chart, Figure 
(85 used to study the feasibility and to expand the basic 
EIU Ussues of computerizing NATOPS charts. As previously 
discussed, each sub-chart was evaluated separately and then 
EU ted into one computer program. However, Sub-Chart IV, 
| шперооГге Clearance, was not utilized in the program since 
| EI X:»ect of take-off is seldom critical for the A-6E air- 


| craft. Each sub-chart family of curves was relatively simple 





Ereproduce using the Least Squares method since the highest 


Brder cof polynomial equation approximated was third order. 


— 





| ihe accuracy of each reproduced curve was verified to be 

| Pecater than ninety-nine percent (99$). A comparison of Sub- 
EI Gross Weight, of Figure 3 with that plotted from 

Bethe approximated equations, Figure 4, will show that an almost 

| exact reproduction of the NATOPS charts can be made by Least 

а р approximations. The ordinate system of Figure 4 is 
meansposed relative to that in Figure 3, since the output 
Martable of Sub-Chart 1 is the input variable to Sub-Chart ll. 

ШЕШЕ the individual curves of each respective sub-chart 

could be reproduced with a high degree of accuracy, high 
errors were generated when approximating coefficients for the 
Mamily of curves in Sub-Charts II and IIT. These errors were 


a direct result of the curves in Sub-Chart II becoming more 
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mom-tinear as the gross weight increased, and of the large 
E Echanges in Sub-Chart III curves as the take-off dis- 
Knees increased. The errors in these sub-charts were 
reduced to less than two percent (2%) by overlapping smaller 
sub-families of curves of more nearly the same linearity and 
воре the gross weights were reduced to an overlapping 
шата фу Of Curves reflecting gross weight changes of ten 


thousand (10,000) pounds, while take-off distance was over- 


ЖЕГЕ ас the mid-curve of the family. The specific overlap 
! points were determined by trial and error until an acceptable 


Error tolerance was achieved. 


Ina computer program which culminated the respective 


' sub-chart programs into one A-6E NATOPS take-off performance 


| computer program is included in Appendix B. Sample program 


EM scare in Appendix E, while Table III shows the relative 


accuracy between computer and manual chart determinations of 


EN Lost distances for representative input parameters. 


24 





TABLE III 


Ae ORF DISTANCE ACCURACY 


x | 
Kas Weight (ibs.) 34,000 . 40,000 | 48 0557000 





| 
Pressure Altitude (ft.) 0 1000 
| 
| 
| 





5,000 0 

Temperature, Deg. (F) 60 100 60 80 

Runway Heading, Deg. SOU S60 360 360 

E mE Direction, Deg. 360 318 042 060 

Wind Velocity 10 20 20 20 

Headwind Velocity 10 ji. AS 10 
EOE Distance: 

x Manual Chart 1,200 2,550 4,400 | 4,400 

Computer 1200 500 4,300 | 4,400 

Accuracy (%) | 100 98 97 100 


нео CRUISE PERFORMANCE 
ШИ 6E NATOPS Cruise Performance Chart - Phase I is 
[| shown in Figure 5, while all four phase charts are included 
in Appendix C. The same basic tachniques used in take-off 
MeeertoOrmance were applied to cruise performance with exception 
of the iterative routine required for the overlying drag count 
ВИНО guideline curves of Sub-Chart II. As in take-off perfor- 
fiance, Sub-Chart I, pressure altitude curve family, was rela- 
tively easy to reproduce with individual curve accuracies 


Meritied as greater than ninety-nine percent (99%). 
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PHASE I — CLEAN AIRCRAFT TRANSFER SCALE 


DATE 15 MARCH 1971 FUEL GRADE: JP-S 
DATA BASIS: ESTIMATED FUEL DENSITY: 6.8 LB./GAL. 


RECOMMENDED 
MACH NUMBER FOR 


ж - MAXIMUM RANGE 
= 5 AT DRAG COUNT 200 
me 
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Sub-Chart I 


ADf1-311 


— — —P — 


Ep Cruise Performance, Phase I. 
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However, high errors were generated when approximating coeffi- 
Ж ЕЕ (5 гот the curve family due to the increased non-linearity 
Eu essure altitude increased. These high errors were 
reduced to less than two percent (2$) by overlapping sub- 
SSS oO curves from zero (J) to thirty thousand (30,000) 
feet and thirty (30) to forty thousand (40,000) feet. 

Im order to achieve accuracies of greater than ninety-nine 
@emeene (995) for the individual curves of the drag count 
Iun sSub-Chart II, a high number of empirical data points 
were used and the approximated coefficients were limited to 
order. This was required to smooth the curves out and 
Bent OSCillations between points. In addition, Sonne points 
in the zero (f) drag count curve were adjusted slightly to 
Maadi tional curve smoothing and better curve reproduction. 
Mile stadjustments were necessary to achieve individual curve 
accuracy, none was required to approximate coefficients for 
ee entire family of drag curves, which resulted in greater 
than ninety-eight percent (98%) accuracy. 

Due to the extreme non-linearities of the guideline family 
ШОЛУ С аг II, the individual curve approximations were 
limited to the range that overlaid only the drag count curves. 
NUES restriction was required to achieve greater than ninety- 
nine percent (99$) individual curve accuracy, and to limit 
the overlapping sub-family of curves to two. The sub-families 
overlapped at the fourth guideline, as counted from bottom to 


ШЕВ іп Figure 4. 
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СС the Hewlett-Packard 9830 Computer Plotting Routine 
mequires that the X and Y axis increase from the ordinate, 
Mm@emeranster scale of Sub-Chart II, Figure 4, had to be in- 
EN w:o reflect values from zero to thirty, vice thirty to 
Exo ihis inversion of the Y-axis caused the approximated 


Ru ume curves to be plotted inverted, Figure 2. However, 


his did not affect the output, transfer scale, since the 
EE "Spoint of entry into the guideline curve family is 
ШЕ poudent upon the initial transfer scale, which is de- 
са rom the zero drag count equation regardless of curve 
orientation. Also, as discussed previously, the iterative 


meOutine required that the approximated drag count family of 


E 5. transposed, Figure 2,so that Best Cruise Mach could 
be determined by comparing input Mach to the guidelines and 


output Mach from the drag counts. An attempt was made to 


keep the drag count curves oriented the same as the guideline 


Ew but the program would not step out of its iterative 


E ute when the input Mach would not pass through the desired 


EN uut curve when intercepting the zero drag count curve. 


For example, as depicted in Figure 4, an input Mach of 0.73 


idoes not pass through a desired drag count of fifty (50). 


Mere nal computer program which culminated the respective 
Bu chart programs into one A-6E NATOPS Cruise Performance Phase 
mecomputer program is included in Appendix D. Sample Program 
membres are in Appendix E, while Table IV shows the relative 


accuracy between computer and manual chart determinations of 


Meranster scale and best cruise Mach for representative input 


parameters. 





ВЕ IU 


Po oe I CIT SET PERFORMANCE ACCURACY 


E odi vs 4 
| | | 
| | Pressure Altitude (ft.) . 33,000 ! 25,000 ЕГЕ ЕЕЕ ООО 
| | | 
Gross Weight (lbs.) | 48,000 | 50,000 | 58,000 | 42,000 
Drag Count 20 | 100 | 150 20 | 





Manual 


| 
| 
Best Cruise Mach 





Computer А .623 .557 


Manual | 

| | о 9.8 (586 14.4 
| 
| Computer oS ES 8 14.5 | 
| | | 
| Accuracy (°) | 99.2 99 98.7 | 99.3 | 
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ТУ. CONCLUSIONS AND RECOMMENDATIONS 


Eun the scope of this thesis project, computerization 
of aircraft Naval Air Training and Operating Procedures 


Standardization Program (NATOPS) flight performance charts 


has been shown to be feasible with a high degree of accuracy. 


The computer programs developed for the A-6E aircraft are 


EE bie tO hand-held calculators, desk calculators, or the 


AS tine aircraft systems computer. The feasibility, pro- 


Mecdures, and techniques discussed are applicabile to any air- 


EET tixed or rotary wing. Utilization of computerized 


MmYALTOPS performance charts for before-flight and during-flight 


mission EEGnugsng will! Significantly reduce the human and 


Manic! errors inherent to visual interpretation and inter- 


BepOtation of the charts. 


ee eduction in chart errors, when coupled with aircrews 
Eme computer access to make airborne aircraft performance 
топ based on all performance variables, rather than on 
Mast experience only, will greatly enhance safety of flight 


in completing any or all mission phases. By enhancing safety 


MO flight through NATOPS performance chart computerization 


meen 1s accessible to all flight crews, the operational readi- 

ness of all Navy/Marine Corps fleet units will be increased. 
Development of NATOPS performance charts computer programs 

Should be expanded to encompass all phases of mission planning 


to include shipboard Bingo (minimum fuel) profiles and tactical 
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mans delivery. These computer programs could be incor- 
ВАО Into hand-held calculators with dedicated read only 
SS ШОМ) for use in ready rooms and in aircraft without 
K systems computers, such as most helicopters and light 
Bu IXed wing aircraft; into desk top computers, capable 
EN потіє cassette tapes to program performance equations 
for all operational Navy/Marine Corps aircraft for use in land 
Ed bp base Operations; and into existing on-board system 
рвет, for such aircraft which have the capacity for 
Computer software expansion. 

Manual inputs to performance equations incorporated within 
Ku systems computers could be significantly reduced by 


винен signals from air data computers for temperature, 


| pressure Рене, and velocity inputs; fuel quantity, sensors 


Во weight inputs; inertial platforms for angle inputs 


ES mangle of attack, climb angles and glide paths; and 


ПЕСН. selection tapes for drag count inputs. 


Eure thesis research should address the interface of 


Mextsting aircraft sub-system output signals as input signals 


| to NATOPS performance chart equations incorporated into on- 


шака ircraft systems computers. 
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AREENDIX N 


mee) SQUARES METHOD OF POLYNOMIAL APPROXIMATION 


The Least Squares method of polynomial approximation is 
Eel in detail in reference [3]. Since data points are 
EN eto experimental errors which may be relatively large 
EN  umuudse, the best curve fit through these points is one 
Euenimizes the errors. The Least Squares method of 


Em mMatioh minimizes the maximum errors between the data 


Ш біп апа the curve by making the sum of the squares of the 
ШЕР minimum. The errors are measured by the vertical 


ЕР ?опса! distances from the points to the line as 


Ed un Figure 6. 





X 


tr QO oh GOn Measurement. 


ШТ Gmm ОЕ the squares of the errors is used since the 
mimetion has no derivative at the origin and this could be a 


miele large error. The Lezst Squares method is also in accord 
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Mene maximum likelihood principle of statistics. If the 
| measurement errors have a normal distribution and if the 
n Ir deviation is constant for all the data, the sum of 
| the squares can be shown to have values of siope and intercept 
| which have a maximum likelihood of occurrence. 

seine development of an an degree polynomial which will 
Mepresemt the best curve fit through empirical data, n+l or 


i reater pairs of data points are required. Let N represent 


EN uUmber Of (х,у) data pairs. The functional relationship 


HESS ta X Fa Xx + .... жа, х (1) 


(2) 
e Y: is used to represent the empirical value correspond- 
¡ing to х; Осе data pair. The Least Squares criterion 


Mequires that 





N 

= 2 2 zo 2 

5 ud 1 F %2 + ° е е е еу x E — 
J | 
2 
N n 5 
P _ й _ DNE - ax. ) (5) 
= E (Y; a aX; а.х. -. Di 


Ве a minimum. 

NI EuSunimum is reached by proper choice of the coefficients, 
ЕА through ал. Therefore, the coefficients are the variables 
ME the problem. At the minimum for S, the partial derivatives 
Sees with respect to the coefficients are zero. Writing the 


lequations for the n*1l derivative equation give: 





N 
/ = = Е - Е К а 
95/92, 0 E 2 (Y; а, а1х. ах; J Ци 
1=] 
М 
35/984 EU a - аух.- po c 
5 = 
52 
а п п 
35/9a. = 0 = f заў “ах; -....- ах: )С- xi ). 


meee each of these equations by -2 and rearranging gives 


EE normal equations to be solved simultaneously: 


2 ne 
a N + ау > х; + а,%х. | an? Xi Dc 
DE c + a ج‎ +a ENS pod cg yx atl = m Ы 
O 1 J 1 Z Í n 1 Qe Due 
2 3 | 4 | mno 2 
aot Xi т ау > х. т аз; 2 х. MU a, 2X; = LX. Yis 
2 (5 
а "E ка am + a r 0. a y x on =s y 
O 1 ] j D Í n 1 QUEE: 


Sere the summations run from 1 to N. 

Smet iacal methods are used to solve these equations simul- 
Laneously and are discussed in detail in reference [3]. These 
particular equations have an added difficulty in that they have 
| the undesirable property known as "ill-conditioning.'" The 
E rs that round-off errors in solving them cause unusually 
marge errors in the solutions, which are the desired values of 
ELE oetricients, a.'s. Up to n=4 or 5 the problem is not too 


1 


Meat, but beyond this point special orthogonal polynomial 
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methods are needed, which is a separate topic in reference [3]. 
Шомемег, from the experimental point of view, functions more 
complex than fourth-degree are rarely needed, and when they 
Е Problem can often be handled by fitting a series of 
polynomials to subsets of data. 

ENUUUStrate the use of equation (5), a quadratic is fit 


Bothe data of Table V. 


TABLE V 


CURVE FIT DATA 


Data Pairs 





Summations for N=6 


2.144 1.9234 3:62.39 ISS 0.74421 


ШОШО л are actually a perturbation of the relationship 


eae - x + 0.2 a (6) 





Maowecurve of equation (6) and the perturbated data are 


Shown in Figure 7. 
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0.6 


0.4 


0.2 


шае 7. Curve fit data plot. 


The summations tabulated in Table V are substituted into 


equations (5) to give the following set of equations: 


6 a en. а 1 Ше 2-55 а) = 9.298, 
5r 2 а 207.5 aj DES a, = в. (7) 
2.5 + DOCU а, + 1.9234 a, = Der ор 


Using numerical methods of solution to simultaneous 


equations yield: 


a, = 0.9986, = -1.0060, a, = 0.2103. 


ш 
Therefore, the equation of the Least Squares approximated 


Polynomial through the data of Table V is: 


y = 0.999 - 1.006 x + 0.210 I (3) 
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iitewcoetricient of equation (8) is slightly different from 
that of equation (6) because of the errors in the empirical 


Gata. 
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TAKE-OFF PERFORMANCE COMPUTER PROGRAM 
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PHASE I — CLEAN AIRCRAFT TRANSFER SCALE 


DATE 15 MARCH 1971 FUEL GRADE: JP-5 
DATA BASIS: ESTIMATED FUEL DENSITY: 6.8 LB./GAL. 
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FUEL GRADE: JP-5 
FUEL DENSITY: 6.8 LB./GAL. 
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CRUISE PERFORMANCE 
PHASE Il] — POUNDS OF FUEL PER NAUTICAL MILE 


DATE 15 MARCH 1971 FUEL GRADE; JP-5 
DATA BASIS: ESTIMATED FUEL DENSITY: 6.8 18./GAL. 
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PHASE iV — FUEL FLOW 


DATE 15 MARCH 1971 FUEL GRADE: JP-5 
DATA BASIS: ESTIMATED FUEL DENSITY: 6.8 L8./GAL. 
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